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Abstract. The colliding wind binary (CWB) systems 7/ Carinae and WR140 provide 
unique laboratories for X-ray astrophysics. Their wind- wind collisions produce hard 
X-rays that have been monitored extensively by several X-ray telescopes, including 
RXTE. To interpret these X-ray light curves and spectra, we apply 3D hydrodynamic 
simulations of the wind- wind collision using smoothed particle hydrodynamics (SPH), 
with the recent improvements of radiative cooling and the acceleration of the stellar 
winds according to a law. For both systems, the 2-10 keV RXTE light curves are 
well-reproduced in absolute units for most phases, but the light curve dips associated 
with the periastron passages are not well matched. In WR140, the dip is too weak, and 
in ri Carinae, the large difference in wind speeds of the two stars leads to a hot, post- 
periastron bubble that produces excess emission toward the end of the X-ray minimum. 



1. SPH Models & X-ray Radiative Transfer 



To model the long-period, highly eccentric colliding wind binary (CWB) systems vj 
Carinae and WR140, we present 3D smoothed particle hydrodynamics (SPH) simu- 
lations of their wind- wind interaction (see Fig. 1), using the same general parame- 
ters a s [Russell et al.|(|2011b|). As improvements upon the previous work of [OkazaH 
et aL| ( [2008| ) and [Russell et al.| ( |2011a|b[ ), the SPH particles are accelerated (using an 
anti-gravity approach) to mimic ayS =1 velocity law, and radiative cooling is explicitly 



treated with the exact integration scheme of jTownsend d2009| ) . 

The SPH visualization program splash ( Price|2007 ) has been modified to calculate 
the absolute flux of the model X-rays from ray-tracing of the emission j = nennAiE, T) 
in the SPH simulations (see Fig. 1), with the emissivity A(E, T) computed from the 



APEC plasma code ([Smith et al.[20 01). Separate energy-dependent opacities, computed 
from windtabs ( [Leutenegger et al. 2010j), are used for the O star and WC star in WR140. 
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Figure 1 . Density (left) and temperature (center) in the orbital plane of the SPH 
simulations, and 3 keV X-rays as observed from Earth (right). Each tick mark is \a. 
The X-ray emission is localized around the apex of the wind-wind interaction region, 
except when the post-periastron bubble occurs in r] Carinae. 



2. Results: Light Curves, Spectra & Hardness Ratios 



Figure 2 co mpares the 2-10 keV RX TE light curve of 77 Carinae ( [Corcoran et al.|2010] 



10]) 



and WR140 (Corcoran et al. 2011) with our model results computed from splas 
Both comparisons are in absolute units, and given this, the level of agreement, apart 
from periastron, is quite remarkable, though there is some cycle-to-cycle variability not 
accounted for in the models. Moreover, figure 3 shows that, for WR140, the model 
and observed apastron spectra also agree, while the model hardness ratio has a similar 
phase variation to that observed, but with an overall softer level that is likely due to the 
limited SPH resolution of the shock-heated material. 

In contrast, near the periastron X-ray minima, both models overpredict the X-ray 
flux. In 77 Carinae, much of this excess results from the fast secondary wind catching 
up with the slow primary wind inducing a shock that blows a hot, extended. X-ray 
emitting, post-periastron bubble into the primary wind. (Curiously, [Parkin et al.| ( |20ir] ) 
find that this bubble does not produce much excess emission in X-rays.) In WR140, 
the nearly equal wind speeds mean there is no such bubble, and so the more modest 
disagreement is just in the depth of the absorption minimum. 

Overall, the WR140 model seems to need only minor tweaking, but in 77 Carinae, 
there are some unknown, fundamentally important effects that are still not understood. 
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^For technical reasons, the RXTE light curve is in flux for r\ Carinae and counts for WR140. We take this 
into account by folding the WR140 model light curve through the RXTE response matrix. 



X-ray Modeling from SPH Simulations 



3 



.X 10 



10 2-10 keV X-ray Light Curve of r| Carinae 



2-10 keV X-ray Light Curve of WR140 




0.5 1 

phase 



1.5 



0.8 



1 1.2 1.4 1.6 l.i 
phase 



2.2 



Figure 2. Model RXTE light curves plotted against the data in absolute units. 
To account for th e Homunculus' absorpt ion in 7/ Carinae, an absorbing column of 
Nh = 5e22 cm"^ ( Hamaguchi et al. 2QQ7| was added. 
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Figure 3. Model RXTE apastron spectra and hardness ratio of WR140, where the 
hardness ratio is the 7.5 keV flux minus the 3 keV flux divided by their sum. 
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